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Abstract 

The method used in [GHS1 to prove that most moduli spaces of K3 
surfaces are of general type leads to a combinatorial problem about 
the possible number of roots orthogonal to a vector of given length 
in Eg. A similar problem arises for £7 in [GHS2]. Both cases were 
solved partly by computer methods. We use an improved computation 
and find one further case, omitted from [GHSlj : the moduli space J-2d 
of K3 surfaces with polarisation of degree 2d is also of general type 
for d — 52. We also apply this method to some related problems. In 
Appendix A, V. Gritsenko shows how to arrive at the case d = 52 and 
some others directly. 

Many moduli spaces in algebraic geometry can be described as locally sym- 
metric varieties, i.e. quotients of a Hermitian symmetric domain T> by an 
arithmetic group T. One method of understanding the birational geometry 
of such quotients is to use modular forms for T to give information about 
differential forms on T\T>. In [GHS1] this method was used to prove that 
the moduli space Tid °f polarised K3 surfaces of degree 2d is of general type 
in all but a few cases. The method works if there exists a modular form 
of sufficiently low weight with sufficiently large divisor. In [G~HSlj, and 
again in [GHS2] where a similar method was applied to certain moduli of 
polarised hyperkahler manifolds, the required modular form is constructed 
by quasi-pullback of the Borcherds form <I>i2. 

A suitable quasi-pullback exists if a combinatorial condition is satisfied: 
there should exist a vector I in the root lattice Eg (or Ej in the hyperkahler 
case) of square 2d, orthogonal to very few roots. This is evidently the case 
if d is large, but for small d the search for such an I invites the use of a 
computer. This was done in both [GHSlj and in [GHS2] by a randomised 
search, relying on the large Weyl group to ensure that in practice no cases 
would be missed. 

Here we present an exhaustive search carried out by the first author. 
For the hyperkahler case the exhaustive search confirmed the results of the 
earlier randomised search, but in the K3 case one previously overlooked 



value of d with a suitable vector was found, namely d = 52. In fact it 
turned out that the randomised search had indeed found this value, and the 
omission of the case d = 52 from [GHS1] happened because the output had 
been interpreted incorrectly. 

Nevertheless the following result is true and has not previously appeared 
in the literature. 

Theorem 1 The moduli space .F2.52 °f K3 surfaces with polarisation of de- 
gree 104 is of general type. 

The paper is organised as follows. In Section [T] we explain briefly what 
the combinatorial problem is and how it arises, and give some more general 
combinatorial problems of the same nature. In Section [2] we describe the 
theoretical and computational methods used to solve it, along with some 
other results obtained in the same way. In Appendix A, Valery Gritsenko 
explains explains how the case d = 52 could have been foreseen without the 
help of a computer. Some of the computer code is given in Appendix B. 
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thesis. He would like to thank Gerard van der Geer for his supervision, and 
the University of Amsterdam for the nice environment it provides. The sec- 
ond author would like to thank the Fondazione Bruno Kessler in Trento and 
the Max-Plan ck-Institut fur Mathematik in Bonn for support, and Valery 
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1 Combinatorial problems and moduli 

In this section we first give a list of combinatorial questions and then explain 
the geometry that originally motivated them. First we fix some terminology. 
We say that L is a lattice of signature (a, b) if L = Z a+b and we fix a bilinear 
form ( , ) : L x L — > Z of signature (o, b). If x G L we refer to (x, x) as x 2 
and call it the length of x. If the length of x is 2 then x is called a root. If 
the roots of L generate L as an abelian group then L is called a root lattice. 
A lattice L is unimodular if it is equal to its dual L v = Hom(L, Z) D L. We 
do not assume that L is always unimodular but for simplicity we do assume 
that L is even, i.e. that x 2 is always an even integer. 

E% denotes the unique even unimodular positive-definite lattice of rank 
8, i.e. with signature (8,0): this is the sign convention of |Bouj and is also 
used in [GHSlj . If n G 2Z then (n) is the rank 1 lattice spanned by a 
vector of length n, and U denotes the integral hyperbolic plane Ze + Z/ 
with e 2 = f 2 = and (e, /) = 1. The symbol © denotes the orthogonal 
direct sum of lattices. If A is a lattice and n G Z, then A(n) denotes the 
same lattice with the quadratic form multiplied by n. In particular, Eg(— 1) 
is the negative-definite even unimodular lattice of rank 8. 

me. - GKS 
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1.1 Combinatorial problems 

Let A be a root lattice (usually it will be Eg or E7) and denote by R(A) the 
set of its roots, i.e. i?(A) = {r G A | r 2 = 2}. The combinatorial questions 
arising in [GHSlJ and [GHS2J are special cases of the following. 

Question 1 Given integers p > q > 0, what are the values of d for which 
every vector of length 2d that is orthogonal to at least 2q roots is orthogonal 
to at least 2p roots? 

More generally we may ask about all possibilities. 

Question 2 Given an even natural number 2d, what are the possible num- 
bers of roots orthogonal to a vector of length 2d ? 

If I G A we denote by R^) the system of roots of A orthogonal to I. We 
denote the answer to Question [2] by P(A,d): that is 

P(A, d) := {m G Z I 31 G A I 2 = 2d, = m}. (1) 

Thus P(A, d) is a finite set of even non-negative integers. We call this the 
root type of the non-negative even integer 2d for the lattice A 

There are some immediate restrictions on what the root type can be: 
for example, if A = Eg then the largest m that can occur is 126, when 
= E-j\ but in that case I G (E-j)^ = A\, so d must be a square. 

Especially for A = Eg, the value of mo(d) = min P(E&,d) is of interest as 
it determines the lowest weight of modular form obtained by quasi-pullback 
(see Equation ([2]) below). If mo(d) = then this form will not be a cusp 
form, so the value of m\{d) = minP(Eg,<i)nN is also significant. We should 
also like to know whether this form is unique. So we also have the following 
questions. 

Question 3 For given d and A, how can we compute mo(d)? 

Question 4 For given m, what is the smallest value d(m) of d for which 
m\{d) < m? 

If in Question U] we replace mi by mo, then the case m = asks for the 
length of shortest vectors in the interior of a Weyl chamber: these are the 
Weyl vectors, which are well known. 

If m G P(A,d) there is a further natural refinement. 

Question 5 How many Weyl group orbits of vectors I with I 2 = 2d and 
#R(l ) = m are there? 
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Some values of m are of particular interest for geometric reasons: for in- 
stance, if 14 G P(E$,d) then quasi-pullback of $12 gives a canonical form 
on J-2d (see Section 11.21 below). This leads us to the following variant of 
Question [TJ 

Question 6 For given m and A, what are the values of d such that m G 
P(A,2d)? 

We can compute the answers to some cases of these questions by the methods 
described in Section [2j 

1.2 Moduli 

The following construction describes several moduli spaces in algebraic ge- 
ometry, including the moduli of polarised K3 surfaces. 

Let L be an even lattice of signature (2,n). The Hermitian symmetric 
domain associated with L is T>l, one of the two connected components of 

V L UV L = {[w] GP(L®C) I w 2 = 0, (w,w) > 0}. 

The group O(L) of isometries of L acts on this union and we denote by 
+ (L) the index 2 subgroup preserving T>l. The action is discontinuous, 
with finite stabilisers, so if V is any finite index subgroup of + (L) then 

F L {T) := T\V L 

is a complex analytic space. In fact it is a quasi-projective variety, hav- 
ing a minimal projective compactification, the Baily-Borel compactifica- 
tion ^(r)*, obtained by adding finitely many curves (called 1-dimensional 
cusps) meeting at finitely many points (0-dimensional cusps). It is of- 
ten preferable to work with a toroidal compactification ^(r), which is 
a modification of J-l(T)* depending on some combinatorial choices at the 
0-dimensional cusps. 

A modular form for T of weight k and character x '■ r — > C* is a holomor- 
phic function F on the affine cone T>* L Cl®C such that 

F(tZ) = r k F(Z) Vt G C* and F(gZ) = x{g)F{Z) V 5 G T. 

F is a cusp form if it vanishes at every cusp. For the cases we shall consider 
the only possible characters are 1 and det(g), and the order of vanishing at 
a cusp is an integer: see [GHS3J. 

The aim of |GHS1] is to show that the moduli space T2d of polarised K3 
surfaces of degree 2d is of general type for most values of d G N. Using the 
Torelli theorem for K3 surfaces one can show that 

^2d = ^L 2d {0 (L 2d )), 
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where O (L) is the finite index subgroup of + (L) that acts trivially on the 
discriminant group L v /L and 

L 2d :=2U®2E 8 (-l)®(-2d). 

Modular forms of suitable weight can be interpreted as differential forms on 
the moduli space provided that they have sufficiently large divisor. There- 
fore, to prove that the moduli space is of general type it is enough to give a 
sufficient supply of such modular forms. There are several technical difficul- 
ties here, one of which is the presence of singularities. A sufficient condition, 
however, was given in [GHSlj . 

Theorem 2 Suppose that n > 9 and that there exists a nonzero cusp form 
F a of weight a < n and character x = 1 or x(d) = det(g), vanishing along 
any divisor H C T>l fixed by reflections in T. Then J~l(X) is of general type. 

The form F a is then used to give many forms of high weight with suf- 
ficiently large divisor, of the form F = F^F^ n _ a ^, and these in turn give 
pluricanonical forms on a smooth model of J~l(T). 

To apply this in specific cases such as T 2d one must therefore construct F a . 
The method used in [GHS1] to do this is quasi-pullback of the Borcherds 
form $i2- This construction first appeared in [BKPSJ. The Borcherds form 
itself was constructed in [Bor] by means of a product expansion, whereby 
its divisor is evident. It is a modular form (not a cusp form) of weight 12 
and character det for the group + (7/2,26)- The lattice 1/2,26 of signature 
(2,26) is 2U © N(— 1), where N is any one of the 24 Niemeier lattices, 
positive definite unimodular lattices of rank 24: see |CSj . For our purposes 
the correct choice of N is 3Eg. A choice of a (not necessarily primitive) 
vector I G Eg of length 2d gives an embedding 

L 2d = 2U® 2£ 8 (-l) © (-2d) 1/2,26 = 2U® 3£ 8 (-l) 

which in turn gives an embedding 

Denote the images of these embeddings by L-2d[l] and V[l] respectively. 

If r G L is a root it determines a Heegner divisor T~L* C P*, given by 
the equation (Z, r) = 0. The Borcherds form vanishes (to order 1) along all 
the Heegner divisors for L = //2,26 an d in particular its restriction to T>*[1] 
vanishes, as needed to apply Theorem [2j However, $i2|x>*[t] ma ¥ wen be 
zero, since if r is a root of //2,26 orthogonal to L 2 d[l] then V[l] C H*. 

Instead we take the quasi-pullback, simply dividing by the equation of 
each such H°, noting that %*_ r = H'. We put 

Ri = {r€ i?(// 2 ,2 6 ) I (r,L 2d [l}) = 0} = {r £ R(E 8 ) | (r,l) = 0} 
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and define the quasi-pullback to be 



F[l] 



$12 



(2) 



]l±re/? ; ( r ' Z ) 



V'[l\ 



This is a nonzero modular form, and one can show that it is a cusp form pro- 
vided Ri 7^ 0. It vanishes along all the Heegner divisors fixed by reflections 



The weight, however, goes up by 1 every time we divide, so the weight of 
F[l] is 12 + ^#Ri- We can therefore show that ?2d is of general type if we 
can find an / G Eg of length 2d with 2 < #R t < 2(n - 12) = 14. Moreover, 
if we can find a cusp form of weight precisely n = 19 then, by a result of 
Freitag [Fr], Tid has p g > and in particular is not uniruled. 

This leads us to Question [TJ with q = 1 and p = 7 or p = 8, for A = £7g. 
In [GHS2j . similar considerations about the moduli of some hyperkahler 
manifolds with a certain type of polarisation lead to Question [1] with q = 1 
and p = 6 or p = 7, for A = Ej. 

2 Solving the combinatorial problems 

The specific combinatorial problems encountered in [GHSlJ and [GHS2] can 
be solved in principle by first bounding d. It is clear that for sufficiently 
large d an I will exist orthogonal to a number of roots in the required range: 
indeed, for sufficiently large d we can find I orthogonal to exactly two roots. 
An explicit bound, followed by a finite calculation, will solve the problem. 
Neither is entirely straightforward, though. In [GHS1] a counting argument 
is used to show that an I & Eg with l 2 = 2d, orthogonal to at least two and 
at most 12 roots, exists (and therefore Tid is of general type) unless 



where Ni(2d) is the number of ways of representing 2d by the quadratic 
form L. The inequality ([3]) certainly fails for large d, but to obtain an 
effective bound on d one must bound NE 6 (2d) and N£> 6 (2d) from above and 
NE T (2d) from below by explicit functions. This is a non-trivial problem 
in analytic number theory but it can be done, and after some refinements 
it gives a reasonable bound of around d = 150. It would be possible to 
resort to direct computation at that point, but there is no need yet. Some 
integers in that range are excluded from the list of possibly non-general 
type polarisations because the inequality ([3]) (or another similar inequality) 
in fact fails. Others can be excluded by inspection, actually producing a 
vector I by guessing the root system R(lg s ). The root systems used in this 
way in [GHSlj were 4^4i, 2Ai © A2, A3 and A\ © Ai- The root systems 
ZA\ © A2 and 2A2 were not tried: see Appendix A. 



in 0+(L 2d ). 



28N E6 {2d) + 63N D6 {2d) > AN E? (2d) 



(3) 
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In [GHS2] a similar procedure was used, although there is an extra dif- 
ficulty caused by the opposite parity of the rank: working in Ej, one needs 
to estimate Njt(2d) from above for some odd-rank root systems R, and this 
problem is not so well studied as in the even rank case. 

In either case, eventually one is left with a residual list of values of d 
for which the problem has not been settled. In [GHS1J it consists of most 
integers between 15 and 60 (for very small d the moduli space is known to 
be unirational) . The residual problem in the hyperkahler case considered in 
[GHS2j is much smaller. 

Now, if we want to be (reasonably) sure that no cases have been missed, 
we do need a computer. Moreover, the methods we now use to solve this 
problem can also be used to give answers to question such as those posed in 
Section 11.11 

2.1 Algorithms 

We begin by representing E$ in the usual way, as the set of points I = 
£ R 8 such that the Zj are either all integers or all strict half- 
integers (i.e. either ij E Z for all i or 2Zj is an odd integer for all i) and 
^2 h G 2Z, with the standard Euclidean quadratic form on M 8 . 

We need a very rough upper bound on Ne s (2d), because we want to know 
whether NE 8 (2d) is small enough to allow a brute-force search for I E Eg 
with I 2 = 2d having 2 < ^R^) < 12. We can easily find such a bound 
by noting that if I 2 = 2d then each of the 8 components 1% of I must have 
if < 2d, so —y/2d <U< y2d, and must be a half- integer: that gives 

N Ea (2d) < (2[2V2d\ +1) 8 (4) 

For d = 52, this bound is about 8 • 10 12 . 

If we are a bit more precise, and note that the components of I are either 
all integers, or all proper (i.e. non-integer) half-integers, we save a factor 
2 7 , giving a bound of about 5 • 10 10 . This is within reach of a brute- force 
search, but it is still high, especially considering that we have to do some 
substantial work for each candidate (compute the inner product with 240 
different vector^!). 

Thus an exhaustive search of all vectors in E% of length < 60 is not 
computationally impossible but it would be cumbersome and would not 
extend to even slightly larger problems such as other cases of Question [TJ 
The Weyl group W{E 8 ) has order 2 14 • 3 5 • 5 2 • 7 = 696729600 and should be 
used to reduce the size of the problem. There are two approaches to doing 
this. 

2 We can be a lot more efficient than that, and skip most of these inner products, but 
even then we still have to compute dozens of inner products per candidate vector. 
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A. Randomised search. This is what was actually done in [GHSlj and 
[GHS2J. Since the non-existence of a vector I gives no information about 
the moduli space, we are willing to accept a very small probability of failing 
to detect such a vector. We therefore choose a large number of vectors of 
length less than 2 • 61 at random and expect that, as the Weyl group orbits 
are large, every orbit will be represented. 

This approach worked very fast, using only a laptop computer and imme- 
diately available software (Maple). A search of twenty thousand randomly 
chosen vectors found all the pairs (d, #ii(/" L )) in the ranges wanted within 
the first two thousand iterations, in approximately two minutes. That is 
fairly convincing practical evidence that there are no more. Unfortunately 
the output was then mistranscribed, leading to the omission of the case 
d = 52 and the erroneous (but not really misleading) statement in [GHSlJ 
that "an extensive computer search for vectors orthogonal to at least 2 and 
at most 14 roots for other d has not found any". 

It is noteworthy that a similar search in the case A = E-j did find some 
cases not discovered analytically, and for which a constructive method of 
finding I is still not known. In other words, some cases of the main theorem 
of [GHS2] still have only a computer proof, although once I has been found 
it is easy enough to verify its properties by hand. 

It is not so easy to estimate the probability a priori that a Weyl orbit 
might be missed. The Weyl group of R^), which is a subgroup of the 
Weyl group of Eg, obviously stabilises I and has order no more than 24 if 
#R(l ) < 12, but in principle the stabiliser of I in W(Eg) could be much 
larger. In that case the Weyl group orbit would be small and more easily 
missed. In practice the randomised method seems to find all the orbits. 

B. Exhaustive search. The first author organised an exhaustive search, 
exploiting the Weyl group by searching a fundamental domain for the sub- 
group H < W(Eg) generated by permutations of the eight components Zj 
and sign changes of an even number of components. This subgroup H has 
size 2 7 • 8!, so index 135 in W(E^): it gives us most of the symmetries, with 
very little effort. 

We say that I £ Eg is in normal form if its components are all nonneg- 
ative (except possibly the first, 1%) and the squares of the components are 
nondecreasing from low index to high index. By acting with an element 
of H, we can translate any I £ E$ to one in normal form: first permute 
the components, so their squares are in order; then make them all (but l±) 
nonnegative, by changing the sign of every negative component (except h), 
and flipping the sign of l\ once for every such change. 

It is straightforward to enumerate the elements of length 2d in Eg that 
are in normal form. For brevity, we will describe this only for the ones 
having integer components (one can get the ones with proper half-integer 
components in a very similar manner). 
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Step 1. For every index i ^ 1, in descending order, we consider all the 
possible values of If we require Zj to be a non- negative integer such that 

• its square, added to the sum of the squares of the coordinates that 
have been chosen (i.e. the I 2 with j > i), does not exceed 2d (otherwise 
I 2 > 2d, for any further choice of coordinates); and 

• (unless i = 8) it is not greater than Zj + i (otherwise I would not be in 
normal form). 

In other words, we let li take any value s£Z such that 

< s < mm j/ m , ^2d- X^f} ■ ( 5 ) 

Step 2. See if 2d — Y3j=2 * s a perfect square m 2 . If so, let l\ take values 
— m and m; if not, discard this choice of coordinates. 

Step 3. Check whether the I so obtained are in E 8 , i.e. whether £| =1 lj € 2Z. 
Discard any that are not in E%. 

We must then filter these enumerated I £ Eg to find the ones with #i?(Z _L ) 
in the required range (2 < j^R{l L ) < 12 for the case considered in [GHSlj): 
this part of the procedure is exactly the same as for the randomised version. 
Since the roots come in pairs ±r it is enough to take inner products with 
a prepared list of positive roots (120 or them), and of course we can stop 
examining I as soon as we find a seventh pair of roots orthogonal to it. 

The first author implemented this search in a high-level programming 
language (Haskell). Without spending much time optimising, this runs fast 
enough (a second or so on commercial hardware, for each of the low values 
of d we are interested in, namely d < 60). The partial use of the symmetries 
of Eg is crucial, though: to go through all the vectors of given length 2d 
would have taken weeks or months for a single value of d. 

This program discovered the lost case d = 52 and therefore Theorem [TJ 
A variant of it for Ef reconfirmed the results obtained by the randomised 
method in |GHS2j . The code used for the case is given in Appendix B. 

2.2 Further results 

The exhaustive algorithm (B) from Section \2. II can be modified to compute, 
in reasonable time, answers to some of the questions from Section 11.11 for 
small values of the parameters. We investigated Question [2] and and Ques- 
tion [6] for small m and d with A = E-j and A = E% . For A = E% we also 
investigated Question [S] for the particular case m = 14, corresponding to 
canonical forms on J-^d- 

Specifically, we have so far computed the root type -P(A, 2d) for A = E-j 
and A = £"8 and d < 150, and the first part of the root type (whether 
m £ P(A, 2d) for 2 < m < 20, say) for larger d, up to about 300 (further 
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for some values of d). This part of the computation is fairly fast and only 
minor changes to the program are needed. 

A little more work, and more computer time, is needed for Question^ We 
must work now with W(Eg), not with H, and we first compute a transversal 
for W(E$) : H (representatives for each of the 135 left cosets of H) and then 
reduce each of the 135 translates of each I to standard form before comparing 
them. 

The outcome counts the number of ways of obtaining a canonical form 
on by quasi-pullback of <3?i2 • There is no assurance either that the forms 
so obtained are linearly independent or that there are not more canonical 
forms that do not arise this way. The results are nevertheless intriguingly 
unpredictable. There are no such vectors for d < 40. There is such a vector 
for d = 40, and also for d = 42, 43, 48 (two orbits), 49, 51-54, 55 and 56 
(two orbits each), 57 and 59. There is no such vector for d = 60, but for 
61 there are three orbits and thereafter the number of orbits drifts upwards 
irregularly. Without further comment, we tabulate below the number ^14 of 
W(E$) orbits of length 2d vectors in E% orthogonal to exactly 14 roots for 
61 < d < 150. 
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A Appendix: d = 46, 50, 52, 54, 57, by V. Grit- 
senko 

In this appendix we find a vector I 6 E% of square 2d orthogonal to exactly 
12 roots in Eg, where d is as in the title of the appendix. (See [GHSlJ 
and [GHS2] for the general context of this question.) We use below the 
combinatorics of the Dynkin diagram of Eg. We take the Coxeter basis of 
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simple roots in £"8 as in [Bouj : 

OL\ Q! 3 OL\ 



«5 



"6 



a 7 



a 8 



0'2 



where (ex, . . . , eg) is a Euclidean basis in the lattice Z 8 and 



1 



(ei + 



«2 = ei + e2, 



1 



(e 2 + e 3 + e 4 + e 5 + e 6 + e 7 ), 
= e k -\ - e k -2 (3 < k < 8). 



The lattice £7g contains 240 roots. We recall that any root is a sum of simple 
roots with integral coefficients of the same sign. The fundamental weights 
cjj of -Eg form the dual basis in Eg = E$ , so (aj, uij) = 6{j. The formulae 
for the weights are given in |Boul Tabl. VII]. The Cartan matrix of the dual 
basis is 



((0Ji, Uj)) 
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We consider the two following cases when the orthogonal complement of a 
vector I in E% contains exactly 12 roots: R{1~e 8 ) = A2 ® 3A\ or A2 © A^. 
(We note that #R(A{) = 2 and #R{A 2 ) = 6.) 

I: d = 46, 50, 54, 57. There are four possible choices of the subsystem 
A2 © 3A\ inside the Dynkin diagram of Eg according to the choices of simple 



(1 3) 

roots of A2, namely A 2 ' 



(7,8) 



(ai,a 3 ), Ap 4) = (02,0:4), 



A 



(5,6) 



(a 5 ,a 6 ) 



or 



(0:7, ag). If A2 is fixed then the three pairwise orthogonal copies of 



A 

A\ in the Dynkin diagram are defined automatically. 

(5,6) 



First, we consider A^'^ 



Moreover A, 



(5,6) 



2^(5,6) - g root system of the orthogonal complement 
of the vector 1$$ = uj\ + W4 + uij E E$ . In fact, if r = Ylf=i x % a i is a positive 
root (xi > 0) then (r, £5^) = x\ + x^ + x-j = 0. Therefore x\ = X/± = xj = 

and r belongs to A 2 5 ' 6 ^ © 3A^' e \ Using the Cartan matrix ([6]) we obtain 
that 1§q = 2-46. Doing similar calculations with the other three copies of 
A2 given above we find 

k,3 = ^4 + UlQ + U> 8 , Z 2 ,4 = W 3 + UJ 5 + U) 7 , l 7>8 = U)i + W 4 + UlQ 



(a5,as). Then 3A\ 



(a 2 ) © (a 3 > © (as) 
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with ftp = 2-50, l\ A = 2-54 and l\ % = 2-57. 

II: d = 52. We consider the sublattice M = A2 ©^2 = (03, 04) © (a6, 07) in 
Eg,. Then M is the root system of the orthogonal complement of the vector 
Im = + 0J2 + ^5 + w§ with = 2-52. 

V.A. Gritsenko, Universite Lille 1, Laboratoire Paul Painleve, F-59655 Vil- 

leneuve d'Ascq, Cedex, France 

valery . gritsenkoOmath . univ-lillel . f r 



B Appendix: Computer code 

Below is the code used to check the combinatorial problem from [GHSlJ, 
and thus to find Theorem [TJ The programs were written in the functional 
programming language Haskell (http: //www.haskellTorg]). The web page 



http://people.bath.ac.uk/masgks/Rootcounts contains links to further 
code and output. 

{-# LANGUAGE TypeSynonymlnstances , NoImplicitPrelude #-} 
module E8 where 



import qualified Algebra. Ring 

import Control . Applicative ( (<$>) , (<*>) ) 

import qualified Data. Vector as V 

import Data. List (intercalate ,nubBy) 

import qualified Data.MemoCombinators as Memo 

import Data. Ratio 

(Ratio , numerator , denominator , (°/ ) ) 

import qualified Data. Set as Set 

import Data.Typeable (Typeable) 

import Math . Combinatorics . Species 

(ksubsets , set , of Size , enumerate , Set (getSet , Set) , Prod(Prod) ) 

import MyPrelude hiding (numerator , denominator , (°/ ) ) 

import qualified Prelude 

import System. Environment (getArgs) 

import qualified Algebra. Additive 

— Some types and helper functions for dealing with 

— "vectors" (implemented as arrays of rational numbers) . 



type Coordinate 
= Ratio Int 



type Vector 

= V. Vector Coordinate 
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— Inner product . 

inp : : Vector -> Vector -> Coordinate 
inp a b = V.sum (V.zipWith (*) a b) 

half : : Coordinate 
half = 1 7, 2 

— Product of scalar with vector. 

1 : : Coordinate -> Vector -> Vector 
1 = V.map . (*) 

instance Algebra. Additive. C Vector where 
(+) = V.zipWith (+) 
(-) = V.zipWith (-) 
negate = 1 (-1) 

zero = V.fromList [0,0,0,0,0,0,0,0] 

— Some data regarding E_8 

delta : : (Eq a, Algebra. Ring. C b) => a -> a -> b 
delta i j = if i == j then 1 else 

— 'e i' gives the i'th standard basis vector of R_8. 
e : : Int -> Vector 

e i = V.fromList $ map (delta i) [1 . . 8] 

— This is the usual integral basis of the lattice E_8. 
basis : : [Vector] 

basis = 
[ 

1 half $ (e 1 + e 8) - (sum $ map e [2 . . 7]) 
, e 1 + e 2 

] ++ map (\ i -> e (i - 1) - e (i - 2)) [3 . . 8] 

roots : : [Vector] 

roots = d8 ++ xll8 where 

d8 = concatMap ((\ [a,b] -> 

[a + b,a - b,b - a, negate a - b] ) . map e . getSet) $ 
enumerate (ksubsets 2) [1 . . 8] 
xll8 = map (\ (Prod (Set neg) (Set pos)) -> 

1 half $ sum (map (negate . e) neg) + sum (map e pos)) $ 
enumerate ((set 'of Size' even) * set) [1 .. 8] 
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— 'posRoots' contains exactly one of every pair 

— (a, -a) of roots. 
posRoots : : [Vector] 

posRoots = nubBy (\ a b -> a == b | | a == negate b) roots 

— Generate elements 1 of the E_8 lattice with the property 

— that 1~2 = 2 d. We need only one element of each orbit 

— under the action of the Weyl group. In particular, we 

— may assume that all coordinates but one (say, the first) 

— are nonnegative, and that the successive coordinates are 

— nondecreasing. We generate exactly one element of each 

— H-orbit, where H is the subgroup of permutations and even 

— sign changes. 

gen : : Int -> [Vector] 

gen d = genlnt d ++ genHalflnt d 

genlnt : : Int -> [Vector] 

genlnt d = map (V.fromList . map f romlntegral) $ go [] where 

— Given the length of a partial vector, compute the maximal 

— new coordinate which does not increase the length of the 

— vector beyond 2 d. 
maxCoord : : Int -> Int 

maxCoord s = floor (sqrt (f romlntegral $ dD - s) : : Double) 

dD : : Int 
dD = 2 * d 

— We maintain a list of coordinates chosen so far, every 

— one together with the sum of squares of the coordinates 

— up to and including that coordinate. 

— The generated vectors are elements of E_8, because the 

— sum of the squares of their components is even, hence 

— the sum of the components as well, 
go :: [(Int, Int)] -> Int -> [[Int]] 

— We have fixed all eight coordinates, 
go f ixed@((_,sq) : ps) 8 

— The vector has the right length; add the relevant 

— solutions (using 'vary'), and continue searching. 
I sq == dD = vary (map fst fixed) ++ lower ps 7 

— The vector has the wrong length, continue searching. 
I otherwise = lower ps 7 

go fixed n = let 

(m,s) = case fixed of 
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[] -> (maxCoord 0,0) 

(c,s) : _ -> (Prelude .min (maxCoord s) c,s) 

in 

go ((m,s + m " 2) : fixed) (n + 1) 

— Lexicographically decrease the given vector, and continue 

— the generation from there. 

lower :: [(Int,Int)] -> Int -> [[Int]] 
lower [] _ = [] 

lower ((x,s) : ps) n 

I x == = lower ps (n - 1) 

I otherwise =go((x-l,s+l-2*x) : ps) n 

vary : : [Int] -> [[Int]] 
vary (x : xs) = if x == 
then [0 : xs] 

else [x : xs, negate x : xs] 

— For vectors with all coordinates half -integers , we work 

— with the doubles of the coordinates. 
genHalf Int : : Int -> [Vector] 

genHalflnt d = map (V.fromList . map (% 2)) $ go [] where 
maxCoord : : Int -> Int 
maxCoord = Memo . integral m where 

m s = f $ floor (sqrt (f romlntegral $ dE - s) : : Double) 
f k = if odd k then k else k - 1 

dE : : Int 
dE = 8 * d 

go :: [(Int, Int)] -> Int -> [[Int]] 
go f ixed@((_,sq) : ps) 8 

I sq == dE = filter e8 (vary $ map fst fixed) 

++ lower ps 7 
I otherwise = lower ps 7 
go fixed n = let 

(m,s) = case fixed of 

[] -> (maxCoord 0,0) 

(c,s) : _ -> (Prelude. min (maxCoord s) c,s) 

in 

go ((m,s + m " 2) : fixed) (n + 1) 

— Decides whether a given vector is an element of E_8 
e8 : : [Int] -> Bool 
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e8 = (== 0) . flip rem 4 . sum 



lower :: [(Int, Int)] -> Int -> [ [Int]] 
lower [] _ = [] 

lower ((x,s) : ps) n 

I x == 1 = lower ps (n - 1) 

I otherwise =go ((x-2,s+4-4*x) : ps) n 

vary : : [Int] -> [[Int]] 

vary (x : xs) = [x : xs, negate x : xs] 
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